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I. MODEL DEVELOPMENT OF THE CONTINUOUS-DATA IPS 
CONTROL SYSTEMS 


The objective of this chapter is to develop the dynamic equations and 
the mathematical model of the continuous-data IPS control system. The IPS 
model considered includes one flexible body mode and is hardmounted to the 
Orbi ter/Pal let. The model contains equations describing a torque feed-forward 
loop (using accelerometers as inputs) which Will aid in reducing the point- 
ing errors caused by Orbiter disturbances. 

The equations of motion of the IPS system are written as 

MR + DR + KR + Q/Rdt = F ( 1 - 1 ) 


where 



6 . : 

i 

X, Z : 
M : 

D : 

K : 

Q : 

F : 



X 

Z 

n 


: translations of the Orbiter 

attitude of the Orbiter 
attitude of the instrument 
accelerometer measurements 
7 x 7 mass matrix 
7><7 damping matrix 
7 x 7 stiffness matrix 
7 x 7 integral control matrix 
7 x 1 generalized force vector 


0 - 2 ) 


The elements of M are: 
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other elements of M are zero. 


+ r bx C *» 


The elements of D are: 


d 

d 

d 

d 
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43 = 

d w = 

K r - 19 

47 = 

K r°rg 

= -76.4 

55 = 

2c a) 

^X X 

= 44 

66 = 

2 «z“z 

= 44 

77 = 

2 Vb 
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All other elements of 0 are zero. 
The elements of K are: 


% * k w ■ K P = 70 ■ 000 

k i .5 = " m i <r b 2 C ‘*' + r bx S f )u ay = 

2 

k., = m.(-r. Sd> + r, C<b)o> = 
Hb I bz Y bx Y ay 

S? " Vss ' ’ 272 

k 55 = “x = 986 
k 66 = »* = 986 

k 77 = "b = ,32 


-2.23xl0 6 

-3.87x10 6 


All other elements of K are zero. 
The elements of Q. are: 

943 “ 944 = K i * >-' x ' 0 5 

947 = K i°ss “ - k11 

AH other elements of Q are zero. 
The generalized force vector Is 
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(1-3) 


Given the elements of the matrices M, D, K and Q, the equations of 


motion of Eq. (l-l) are written as 





m ll X s + m 13®s + m )lt 0 i “ F x 
m 22 Z s + m 23®s + m 24®i ’ F z 
m 3l X s + m 32 Z s + m 33®s + m 3^®i “ r c.x F z ’ r cz F x 

m <il X s + "42 2 s + "Vs + m W®i + K A + K r®i + VVg" 

+ K p®s + K p®i + %* + k l|6 Z * k ‘.7 n + K i/ 0 s dt 

+ K./e.dt + K.a ss /ndt - 0 

+ d 55 X + k 55 X — ^ 

m 62 Z s + m 63®s + m 66 Z + m 67'' + d 66 Z + k 66 Z = 0 


m__n + d _,n + k_,_,D = h F + h F 
77 77 77 zc z xc x 


(1-4) 
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0-9) 

( 1 - 10 ) 


Without the nonlinear wire-cable and flex-pivot torques, the control 
torque is expressed as 

T c = " K A - Vi - K P 0 s - K p 0 i ‘ K i/ 0 s dt • K i/ 0 i dt 


“ k 45 x " k 46 Z 


0 - 11 ) 


where K p , and K. are the rate, position, and integral constants of the 
controller, respectively. 

The nonlinear torque due to the flex pivot and wire cable can be lumped 
into one operator N which operates on 0.. Thus, the component, -N(0.), 
should be added to Eq. (1-11); i.e., 

T c ■ - K r®s - K r®i - K p®s - K p®i - V e s dt ‘ K i/®i dt 


- k 45 X - V * N(0 i> 


Solving for X g from Eq. (1-4), we have 
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Solving for Z g from Eq. (1^5), we have 

m 23 •• m 24 " 1 

Z = - 0 - — 0. + F — L - 

S m 22 5 m 22 1 2 m 22 

Substitute Eq. (1-13) and (1-14) into Eq. (1-6), we have 

nu.m. m,,m, nu, m.,„rru., m„„m„ , m 

11 e s 


(1-14) 
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(1-15) 


The last equation is rearranged to the following form: 
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(1-16) 


Substitute Eqs. (1-13) and (1-14) into Eq. (1-7), and making use of Eq . (1 
we have 
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Simplifying, the last equation becomes 
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m. . 

+ VV + k 47 n + K i a ss/ ndt ' T c = F x 

Similarly, Eq. (1-9) is written as 
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(1-18) 


+ m 63 0 s + m 66 Z + m 67 n + d 66 Z + k 66 Z = 0 
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Equations (1-16) , (1-18), (1 — 19) > ( 1 —8) and ( 1 — 1 0 ) are now written as 
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( 1 - 20 ) 
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These last five differential equations are rearranged so that a block diagram 
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can be constructed. 
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The control torque T is given by Eq. (1-12). 

c 


(1-37) 


The block diagram which portrays the differential equations of Eqs . 
(1-33) through (1-37) is shown in Fig. 1-1. The nonlinear element which 
represents the nonlinear torque due to the flex pivot and wire cable is also 
included in the block diagram. The signal flow graph representation of 
Fig. 1-1 for the purpose of evaluating the determinant is shown in Fig. 

1-2. The dynamics of q is eliminated since they do not enter the determinant 
of the signal flow graph. 

Applying Mason's gain formula to Fig. 1-2, the determinant of the 
signal flow graph is evaluated as follows: 
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Rearranging both sides of Eq. (1-38), we have 
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Dividing both sides of the last equation by the terms that do not contain 
N, we get the equivalent 1 inear. transfer function that the non 1 i near element 
N sees, 

M s 
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For the system parameters given, G = G ; thus, Eq. (1-42) is simplified to 
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(1-43) 
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Substitution of the system parameters into Eq. ( 1-43) with 

M s = 7056977.95 

M k = M. = 25992.12 

M n = 10074.90 

M = 0.926522 

u 

M p = -4.537982 

M r = 0.042982 

M = 0.02489 
v 

we have 

, v 1 , 002077 x 1 0~ 4 s (s 2 + 44s + 986 ) . 

G eq u; s b + 45.96676 s 4 + 1 1 07.4783 s d + 2257-769s z + 7374.0808s + 10828.48765 

(1-44) 

It is interesting to note that the system shown in Fig. 1-1 is of the 

11th order. However, since the dynamics of q do not enter the picture on 

the stability of the nonlinear system, and since G = 

X 

seen by N is only of the 5th order. 


G z , the transfer function 


II. PREDICTION OF SELF-SUSTAINED OSCILLATIONS OF THE CONTINUOUS- 
DATA IPS SYSTEM WITH WIRE-CABLE AND FLEX-PIVOT NONLINEARITIES 


The equivalent transfer- function G e ^(s) in Eq. (1 _ 44) can be plotted 
in the gain-phase coordinates together with the plot of -1/N of the nonlinear 
element for stability analysis. The describing function of the combined flex- 
pivot and wire-cable nonlinearity using the Dahl model has been derived in 
[1] . In this earlier report a simplified model of the IPS system was obtained 
by assuming that all but motion about two of the seven degrees of freedom of 
axes are negligible. Only motion about the scientific instrument axis and 
the mount rotat ion were cons idered . 

For the purpose of comparison, the equivalent transfer function that 
the nonlinearity sees for the simplified IPS in [1] is repeated as follows: 


G ( s ) " 


eq 


0.001 3946s (s 2 + 0.0012528s +0.0036846) 
+ 16. 7s 4 +222. 6s 3 + 279.806s 2 + 5.857s 


+ 5-13855' 


(2-1) 


when the integral control constant Kj is 10 6 . The zeros of this transfer 

function are s = -0.0006264 + jO. 06069764 and s = -0.0006264 - jO. 060676/4. 

The five-poles are at s = -1.38135, s = -0.0031956 + jO. 135895, s = -0.0031956 

- j0. 135895, s = -7.65618 + jl 1.946, and s = - 7 . 656 1 8 - jl 1.946. It was 

concluded in [1] from the gain-phase plot of G (s) versus the - 1 /N (A) plot 

eq 

that the simplified IPS may have a sustained oscillation that is characterized 
by the frequency of 00 = 0.16 rad/sec. Depending on the values of the para- 
meters of the wire-cable and flex-pivot nonl ineari ties , K^. and H^, the 
amplitude of oscillation of e, which is comparable to 0. in this report, 
lies between 3*10 -8 to 3 X 10 -6 for Kj = I0 6 . 

The equivalent transfer function of Eq. (1-44) is plotted as shown in 
Fig. 2-1. It is noted that the IPS model used in the present work does not 









































have the same system parameters as the simplified IPS system studied in [lj. 
Thus, the natural frequencies of the two systems are different. From Fig. 
2-1 it is observed that the G (s) and -1/N loci would intersect only at 
very low frequencies and when the amplitude of oscillation is extremely 

- 3 

small. For instance, at oo = 10 rad/sec, the amplitude of oscillation 

is approximately 10 10 , and for all practical purposes this can be regarded 

as zero. Thus, the IPS model considered here would not exhibit sustained 

oscillations due to the wire-cable and flex-pivot nonlinearities. 

It is of interest to investigate the poles and zeros of the transfer 

function of Eq. (1-44). The poles and zeros of G (s) of Eq. (1-44) are 

eq 

tabulated as follows: 

zeros: s = 0, -22 + J22. 4053565, -22 - j 22. 4053565 

poles: s = -1.66963, -0.123519 + j2.5232 

-0.123519 - j' 2.5232 
-22.025 + j 23. 0467 
-22.025 - j 23. 0467 

The natural frequency of the dominant poles of the G (s) of Eq. (1-44) 

eq 

is found to be at 2.5232 rad/sec. This is much higher that that of the 
transfer function of Eq. (2-1). Furthermore, the two zeros listed above 
are very close to two of the poles. 
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III. DIGITAL COMPUTER SIMULATION OF THE CONTINUOUS-DATA 
NONLINEAR IPS CONTROL SYSTEM WITH DAHL MODEL AND 
ONE FLEXIBLE BODY MODE 

The IPS control system with one flexible body mode and the nonlinear 
flex pivot torque modelled by the Dahl solid friction model (i = 2) is 
simulated on the digital computer. A block diagram of the IPS system is 
shown in Fig. 1-1. For the Dahl model with i = 2, the following parameters 
are used: 

y = 9.2444 x 10 4 
T fro = 2.25 x lo' 3 N-M 

The objective of the computer simulation of the system is to verify 
and/or clarify the analysis results obtained by the describing function 
method in Chapter II. Although the describing function method analysis 

i * • 

in Chapter II does not yield precise results for possible amplitudes and 
frequencies of sustained oscillations, it does indicate that for the sets 
of system parameters used if sustained oscillations were to exist the 
amplitude of oscillation would be approximately less than lO' 3 radians and 
the frequencies would be less than 10 -3 radians/sec. 

Since the significant time constant of the IPS system is very long, 
it would require a long computer simulation time to verify the existence 
or nonexistence of a limit cycle.’ Normally, the possibility of sustained 
oscillations could be checked by using two computer simulations, one with 
initial conditions that correspond to above and the other with initial 
conditions that correspond to below the value of the predicted amplitude 
of oscillations. If a sustained oscillation exists then the simulation 
with initial conditions below the predicted amplitude will result in a 
response that oscillates near the predicted frequency and increases in 



amplitude. Also, the simulation with initial conditions that correspond to 
an initial amplitude above the predicted amplitude will result in a response 
that decreases in amplitude and oscillates near the predicted frequency. 
However, since the analysis of Chapter II does not predict any specific 
sustained oscillations, the method described above must be changed. 

Since the analysis of Chapter II provides no specific candidates for 
amplitudes and frequencies of limit cycles (or the amplitude and frequencies 
are all extremely small), it may be necessary to check for the absence of 
limit cycles. The absence of limit cycles can be verified by showing that 
the responses of the system due to any initial conditions actually decay to 
zero as time increases indefinitely. 

Again, since the time constant of the IPS system is very long, the 
simulation of the system response by a digital computer may involve excessive 
computer time. Since if sustained oscillations were to exist, the possible 
frequencies of oscillations would be less than 10" 3 rad/sec, the simulation 
time required would be in excess of 2tt x 10 3 sec. Also, a double-precision 
program with 30 digits of accuracy would be required, since the amplitude of 
oscillation must be less than 10 9 radians. Instead, the results from the 
simulation will be used to verify that the osci 1 lations of the system are 
damped out at a mode that is practically independent of the nonlinearity. 

As a means of exciting the system, an initial perturbation of G . ( 0) = 

10 9 radians was applied to the system. Two simulations of the system were 
completed for two extreme sets of wire cable parameters (K^. and . 

The responses of 0. for both of these cases are shown in Figs. 3 _ 1 and 3“2. 
These responses indicate that the influence of the wire cable parameters on 
the system response is very insignificant.- In addition, the frequency of 
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oscillation is approximately 2.5 rad/sec, which is approximately 10** times 

greater than the maximum possible frequency of oscillation. Also, the envelope 

of the response is decaying with a mode of a = -0.1239, where 1/a is considered 

as the time constant. Hence, the system over the 50 sec interval of simulation 

time has a response due to a pair of poles at s = -0.1239 + j2»5 and s = -0.1239 

- j2.5, which is approximately one of the pairs of poles of G (s) . Thus, 

eq 

the system oscillates at the dominant mode of G e ^(s) as expected, if the system 
had no limit cycles and the nonlinearity were replaced by a unity gain element. 
This implies that the nonlinearity has very little influence on the the system 
response under the simulated conditions, and that no limit cycles exist. 

Since the size of the feasible amplitudes predicted is so small, at least 
15 digits are needed to carry out the 50 seconds of computer simulation. 

The listing of the computer simulation program is given in Table 3“1* 


Figure 3-1, Computer simulation program of cont jnuos-data IPS system. 



SD2) 

AK77, HZC t HXC 


1 

C SET 
C a A3 3 


c set 
c set 

C SET 

C SET 
Z SET 

C DBF; 


c sar 



NDi:i=H 

DO 1 1=1 , NDT’1 

IF (I. 30. 3) 30 TO 1 

5f(I)='?.'33 

CONTINUE 

SYSTEM CONSTANTS 

nlTC ?A ivq A .3 « 

a 1 1 =9 . 1 2 74 3 3 1 3 = -2 . 2 6 7? * 3 1 4=2 . 27 E3 

.122 = 3 1 15323=1 . 5654*a24=3. 92F3 

331=3 133332=1233333=7, 06E6?M34«2. 6734 

.1 4 1:= 3 1 4 £ 34 2 = 1 2 4 53 43 = M? 4 SMU 4 = 1 . 0 3 E4 

151=- 1. 7,0 3 15 3 = 9.297-1 0 M c 5= 1 . 70S M 57= - 1 . 1 3 E- 3 

352 = - 1 . 30.S363 = -4.72ETf,M66 = M55$M67=1 . 37E-3 

377=355 

DAMPING CONSTANTS 

D 4 3 = r . 9 7 E 4 3 D 4 4 = D 4 3 S D 4 7 = - 7 6 . 4 EC 

055=4 4. 305O66 = D55$D77=, M5E0 . 

STIFFNESS C 0 M 3 T A N T S 

AK4 3 = 7 . 74 5 A K 4 4= AK 43$ AK.4 5 = -2 . 23 E6 

A K4 6 = - 3 . 37 3 6 5 A K 4 7=- 2 . 7 2 E 2$ A K 55= 9 . 86 E 2 

A K 6 6 = 9 » 3o72SAK77= 1, 32E2 

P ID CONST A M S 

A KR= D 4 4 $ 31 G HG = 047 /A K 3 

,zi.K P = A K44 .5 S I GS S - A K47/A K p 

043=1 . 1 35 oAKI=043$Q44=Q43 

Q47=- 4. 27 33 

0 IS TA ’IC 73 FROM CENTER OF MASSES S CENTER OF GRAVITIES 

TO ONE IF THE CREW MOTION IS NOT INCLUDED IN THE MODEL 

dZC = 1 . E 9 $ H X C = 1 . E9 

RCX = 1 . E 0 $ 2 C Z = 1 . EO 

N-7 AGRFGAT 3 IASS ELEMENTS 

3I = M34- 331* Ml 4/M*» 1- K?2*”2 U/h22 

3 K = H 4 3 - M 1 3 * M 4 1 /M 1 1- M 4 2 * M 23/M 22 

3S = M3 3- 31 3* 33 1/Ml 1- M?2*F23/M22 

3 N = M 4 4- 3 1 4* M4 1/MI 1- X 42* M2 4/M 22 

M U =M 5 3 - ,3 51* 11 3/31 1 

3V = -3 51* 3 1 4/3 1 1 

3 P = - 3 5 2 * '1 2 3 / 1 2 2 + 3 6 3 

1 R = -Mo 2 * 32 4 /3 2 2 

ERROR WEIGHTING FACTORS 

DO 5 1=1 . N 913 

DEFY (I) =1.7 3/11 .39 

E KGS FONTS 07 PARAMETERS 

P RilT ( 1 } =9. EO 

PRUT (2) =0.5 S3 

PRtfT 3) =1 . E-2 

PRtfT ( 4) =1. 7-1 4 

PRM T (5) =0. .7 9 

PRMT (6) =2. 53-1 

P RM T ( 7 ) = ? R 1 T ( 1 ) 

CALL RKS5 (PR3T # Y, DER Y , NDI M , I HLF # FCTY ,0UT?,7.Tm 

STOP 

END - 



Table 3” I (continued). 
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SUBROUTINE DU IP (T IMF , Y , PSRY , IHLF, NDTfi. PRMT) 
DIMENSION Y [ndI M) # Y ( NDIM) , PRHT (7) ' J 


— -- - — — - - — \ - — * i w •<*» * i lx i « 4 t.ui a 

DJ3L=TIflS«P3HP (71 .. 

IF f A3 S TIME-0. SO) . L F. 1.E-1 0) 00 TO 
I F j A3 S (°RMT ( >) -DEI) .. LE. 1. E-10) GO 
RETURN 

1 ?RM! (7) = I-:‘1S 

WRITE (6, 10) IHL? ,T IME ,Y (3) .Y (8) 

10 FORMAT (16, 1PE10.4,1P§ET1 5) 

RETURN ' 

. END • 


1 

TO 


T 


SUBROUTINE FCTY (TIME, STA TS, STDOT, KEEP, N DIM ) 
IMPLICIT REAL ( M ) 



HZC 


Z THIS SUBROUTINE ONLY UPDATES THE DERIVATIVE OF THE STA*E VECTOR. 
C T HE UPDATE DEPENDS ON THE INPUTS, TIME, STATE, & K^EP." 

DIMENSION STATE (NDIM) , STDOT (NDIM) 

?X (TIME) =0. ED 
FZ (TIME) =). E 0 

C 10 FORMAT (16, 1 PE 10. 4,1 P2E1 1.3) 

C 'lODiL THE WIRE NONLINEARITY 
3GNED7=1.S0 

IF (STATE (3) .LT. 0. E A ) SGM EDT=-1 . ED 
TWC=S GNEDT^HUT* AK HASTATE (3) 

Z MODEL FOR THE DAHL NONLINEARITY (1=2) 

I F (Jj'T IME. ST. 0 . EO) . AN D. (KEEP. EO. 0) ) GO TO 1 


30 TO 
11=1 
R 1 = R 0 - S . 


I he or 


TFP= jfSGNEDT>Rl/ (1 .E0-G1 
IF (IFP, GT. T F ? 0 ) TFP^^pn 
IF (TFP.LT, -X?P0) TFP=-TFP0 
IF (II . EQ. 0) GO TO 3 
R0=T?P/TFPO 
1 1=0 

3LAS!=S TATE (3) 

•30 T’O o 
TPP=3 . FD*TF? 

ENN = T PP+ TNT 


STDOT 

STDOT 

STDOT 


* (STATE (3) -SL AST) *R1 ) )*TFPC 


) = 3 T A T 3 

[6) 

= ST AT E 

7 

=57 A TE 

3 

) =S TATE 

9) 

= ST ATE 

'10 

1) = STATE (2 


u i (jTii rvr-’-o . uu : \ t i j ■* i\ ^ a ris 14 +Af.4 A r 

1 5) + AKR* (STATE (7) ♦ STA^fR) +S r G?G*ST A^E (6) ) 

STDOT ( 6 ) = - D 77/M l n ESTATE ( 6 ) -AK77/M77TSTATE (1 ) 
l +HSC/DO 7*F3 (TIME) +FX (TINE) *H>:C/D' 7 7 
STDOT (7) =1 , EO/ ( 1 . EQ-MJ*MK/M3AMf<) * (EI/MS/MU*TCM+ (Kl*!441 /(MS*KN 
'Mill- (RCS + M31 /Ml 1 ) /MS) *FX (TIME) + (M1*M42/ (MS*FN** 1 1) + (*'CX 
-M32/M2?) /MS) *?Z (TIME ) . ; 1 



RETURN 

END 



c 

c 


Table 3-1 (continued), 

S OB P.O IJTI N 3 3KS S (P RMT , Y , DEFY , NDIK , IHLF, PCT , 0 UTP , AIJ X) 
DIMENSION Y (JDI M) ,D2?Y (NDIKV ,AUX (8,FDIM) ,A (4) ,B (4) 
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C 

c 


c 

c 


c 

c 

c 


c 

c 

c 


DIMENSION Y(NI 
1 ,C (4) ,PBMT (7) 
DO 1 1=1 , NOIN 

I . . • # A — * ^ > 


c 

c 


1 A OX (3,1) =.0 6666667* DEFY (I) 

X=PR.1T (!) 

XSND= PENT (1+1) 

H=PRMT (1+2) 

PRMT (1+4)='). 

CALL FCT (X,Y,DERY,KESP,NDIM) 

ERROR TEST ■ 

IF (H* (XSND-X) ) 3 3, 37,2 

PRSPAR A TICN S FOR RUNGS- KOTTA METHOD 

;i)=.5 

2) = . 2923932 
3 =1.707107 
4)=«1656657 
1 ) = 2 . 

2 ) = 1 . 

3=1. 

4=2. 

1 =. r ' 

C (2) = . 20 239 32 
0 3) = 1 . 7071 07 
C { 4 ) = . 5 

PREPARATIONS OF FIRST FUNGE-KUTTA STEP 
DO 3 1=1, N D I M 
AUX (1 ,1) = Y (I) 

AUX 2,I)=D*EY (I) 

AUX (3,1 =0. 

3 AUX (6 , 1) =0 • 

I REC= 0 
H=H + H 
I HLF= -1 
I STSP=0 

I EN D= 0 ' 


START 0 s * A RUNGS- KU T T A STEP 
IF ( (X+H-XSNO) *H) 7,6,5 

5 H=XEMD-X 

6 I SN D= 1 

RECORDING 0? INITIAL VALUES OF THIS STEP 

7 CALL OUT? (X,Y , D ERY , I RFC , NDIM , PRMT) 

IF (PR MX (1+4) ) 40, 3,4 0 

3 ITS3T=0 
9 I S TE ? = I S X3? + 1 


START OF INNERMOST FUN3F-KUTTA LOOP 
J = 1 

10 A J = A (J) 

3 J = B ( J) 

.**■ J ~ P l J I 

DO 11 I=1,NDIM 
R 1 = H * DE ?. Y (I) 

R 2 = A J * ( R 1 - *3 J* AU X (6,1)) 

' Y(I)=Y I)+R2 
R 2 = R 2 + R 2 + R 2 

AUX (6,1) =A'I X (6,1) +R2-CJ*R1 
I? (J- 4) 12,15,15 
J=J + 1 

IF ( J- 3) 13, 14, 13 '* 

X=X+.5*H 
14 CALL FCT (X,Y,DERY , KEEP, NOIM) 

-GOTO 10 

END OF INNERMOST P.UNGS-KUTTA LOOP 


RISPEODUC^'L-Tr OP l 


11 

12 

13 


Table 3-1 (continued). 
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C 

C 

c 

c 


15 


15 

17 


13 


TEST OF ACCURACY 
IF (ITSST) 16,16,20 


IN CASE ITEST=0 THSFE IS NO POSSIBILITY FOP TESTING OF ACCURACY 
uj 17 1=1 f * Dj. H 

AUXJ4,p*Y(I) 

ITEoT=1 


I STEP -IS TCP 1 1 STEP-2 
I HL F= I HLF+1 
X=X-H 
H=. 5*H 

DO 10 1=1 /'HIM 
Y (I ) = AUX (1,1) 


DERY (I ) = AUX (2.1) 
19 AUX (5.1) =AU>: (3,1 
GOTO 9 


) 


C 

c 


IN CASE ITEST=1 TESTING OF ACCURACY IS POSSIBLE 
IMG D= IS TEP/2 

IF (ISTEP-IMOD-IMOD) 21,23,21 
CALL ? C T f.C , Y , D E 3 Y , K 2 F P , V D I M ) 

DO 22 1=1 , NDIM . 

AUX (5,1) = Y(I) 

AUX (7 ,T =DSRI (I) 

GOTO 9 

COMPUTATION OF TEST VALUE BELT 
DELT= 0. 

DO 24 1= 1 , N DIM 

D EL T= D EL T + A UX (S .1) *A3S (AUX (4 ,1) -Y(I) ) 

IF (D3LT- PRMT (1 + 3) ) 2 8,23, 25 . 

ERROR IS 'TOO GREAT 
IF (IHLF-10) 26,3 6, 36 
DO 27 1=1, NOIii 
AUX (4,1) = AUX?5 , 1) 

ISTE?=ISTS?+ISTSP-U 
X = X - H 
I END=0 
GOTO 13 

: RESULT VALUES ARE GOOD 

23 KBS* 5 =0 

CALL" FCT (X, Y , DERY, KEEP, N DIM) 

DO 29 1=1 , N D I -A ■ 

AUX (1 ,1) sYjfl) 

AUX (2,1 =OERY (I) 

AUX (3 ,1 =A'J X (6,1) 

Y (I) = AUX (5, I) 

29 DEFY (I) =AUX (7 , 1 


20 

21 

22 

23 

24 


25 

25 

27 


C 

r* 


34 


30 

31 


32 

33 
35 


CALL OUTP (X-tl , Y, OS? Y , IHLF , NDIM , PRMT) 
IF (PENT (1 + 4) ) 40 ,30, 4 > ’ 

DO 31 1=1, NDIM 
Y m = AUX (1,1) 

DERY (I) = A'JX (2 , 1) 

TDT?r;-2HLF 


I? 

15 


? (KEEP . 3 0. 2) GO TO 39 
F (TEND) 32,32, 39 


INCREMENT GETS DOUBLED 
I HI. ?= IHLF- 1 
I STE ? =1 STEP /2 
H = II + H 

IF (IHLF) 4, 33, 33 
IHOD= ISTEP/2 
IF 
I; 

I HL ! 

ISTE?=I3TEP/2 
H = H + H 
GOTO 4 


:F (IS TH? -I ’100- 1 MOD) 4.34,4 
:F ( D3 LT- • 0 2 * ? RM T (1 + 3') ) 35,35, 
:hlf=ihl?-i . 
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Table 3”1 (continued), 


•24 


C 

c 


RETURNS TD CALLING PROGRAM 
35 ‘ IHLF= 1 1 

CALL ^CT (X, Y, DSRY ,KEEP,NDIM) 

37 imp* 12 
GOTO 39 
33 I HLF - 1 3 

33 CALL DUTP(X,Y ,D3RY,II!L P ,NDIK,PRMT) 
40 3 3 T URN . ' ' 

2ND • ■ ’ 


PRO SR VI PLOTS F ( INPUT /TAP35=INPTJ T , GPH2 , TAPF1 =G?H2, IPSD2, 
1 TAPE6=I?3D2) 

“ ALL GR A P3 ( 1 , "T HETI " , 0. 1 0 9 . 0 , 1 , 2 , "TH DOT'* , 5 , » DUN K 1 " , 1 


1 , " N " , 132) 

END 


SUB ROUTINE GRAPH (!7UHCOL , HEAD, ATM IT , FIN AL , NUMPRT , NIJN P 1 , HE AD 1 , 


NUN ?2 / HE A 02 ,NUUSKP, WATT, WITTY) 
DIM EM3I0 'I D AT A ( Y 0 9) , TIN ES (1 C 00) , 


, r CSTOPE (81), DATA 1 (10 00) , 

* DATA 2 (1001) , CARD (30) , v (15) 

Til 3 GRAPHING PARAMETERS ARE THE FOLLOWING: 

THE NO. THF COL, OF THE DATA FILE TO BE PLOTTED 
(AND PRINTED) 

HEADING FOR THE PLOTTED VARIABLE 
AIST = THE INITIAL PLOTTING TIME 
FINAL = THE FINAL PLOTTING TIME 


NUMCOL- 
HEAD = 


(IN QUOTE MARKS) 


NO. OF OTHER VARIABLES TO EE PRINTED ALONG SIDE 


H0MP3T- T4: 

THE GRAPH (0,1,2) 

NUMP1 = NO. OP THE COL. OF m H2 DATA PILE TO BE THE FIPST 
PRINTED VARIABT E 

HEAD 1 = HEADING F0 D THE FIRST PRINTED VARIABLE 

NUMP2 *= NO. OF THE COL. O* THE DATA FILE TO BE THE SECOND 


HE AD 2 = HEADING FOR THF SECOND PRINTED VAFI ABLE 

SUMS KP= NO. OF LI NFS OF DATA FILE TO BE SKIPPED BETWEEN 

SUCCESSIVE PRINTS (THIS ALLOWS THF GRAPH TO BE COMPRESSED) 
WAIT - "W" CAUSES m HE PROG . TO WAIT FOP A CARRIAGE RETURN INPUT 
BEFORE IT STARTS TO PRINT. THIS ALLOWS TIME SO THE PAPER 
CAN BE POSITIONED TO THE TOP OF A NEW PAGE. 

NUMT TY= THE NO. OF CHARACTERS PER LINE OF , TTY OUTPUT. (80,132) 

I? 30 IS USED THFN NO ADDITIONAL COLS. ARE PRINTED, 

IN PIT format re; 


1 FORMAT REQUIRED: 

FORM AT (6X,1?E10. U, IpIQEll . 3) 

f,Y= frrST.D UTfiTH ^5'! R V 


6 X= THiS FIELD WIDTH CAN BE' USED TO STORE IHLF OP SOKE OTHER FARM. 

1 PEI 0.4= THIS FIELD IS USED TO STORE THE INDEPENDENT VARIABLE (TIME) 
IP 1321 1.3= THESE FIELD WIDTHS CONTAIN UPf TO 10 COLS. OF DATA 
FOR 10 DIFFERENT VARIABLES 

Q **#£*#** **£^*>!e** ❖*$**£•£* *<##******#>M.s **#*>«** *********£***#* * * * * * if. * * * * * * £ 

IF (WAIT.NE. "W") GO TO 1001 
PRINT 1032 


1002 


FO? 132 '(OR 80) ■ CHARACTERS 
'* TYPE, CARRIAGE RETURN ' TO _ BEGIN OUTPUT'*) 


FORMAT ( " IS 


-> n 




1033 

1031 


ID 


2D 

33 


READ (5,1030) (CARD (L) , L-1 , 3C ) 

FORM A T{ 30 A1) 

ISKIP = 3 C 0 09 
TTY HUM =7 9. 3 

IF(NUMT?Y. EQ. 132) GO TO 10 
TTY MUM =35 . 0 
HUM PRT=0 

AN AM E = ,, TIM E » 

MAX MUM =0 
M A X C 0 L = HU U COL 

IF (NUM ?1 . GT. MAX COL) M AXCOL = NU KP1 
IF] HUMP2 . GT . MAX COL M AXCOL=NU MP2 

READ (6,33 ) fT1 TME, (X (L) ,L=1 , KAXCOL) 

IF (FO? (5 ) . MB. 0) GO TO ^0 
FORMAT ( $ X , G 1 0 .4, lOGl 1 . 3) 




25 


100 

110 

120 

130 


40 

50 

60 


70 

75 


77 

90 


92 

80 


Table 3-1 (continued). 

IF ( (1 , OS - 05 + TIM 5) .LT.AIWI?) GO TO 20 
IS.<i?=I3KI?+1 

IF (I3XIP. L2. N T JM5KP) GO TO 20 
I3KIP=7 

IF( (TIHE- 1.0 3-° 5) .GT. FINAL) GO TO 100 

maauum- ma x:iua+i 

TIMES (MAT TJ A) = ? I'M E 

d \? a { a a x u.'i )=:■:{ no a co i, > 

. DATA) (MAT NO ft) (N II M pi ) 

DAT A2 (ilAXNUM) =X (NOMP2) 

GO TO 20 

I? (N0MP3T. 30.0) WRITE (1,110) ANAME, H-EA D 
forma? (" " . a io , 3x , ai oy 

IF (NUMPRT, 30, 1) WRITE (1 . 12C)A NAHE ,HEA D,HEAD1 
FORMA? ( " ",A1C,3X,A1P.33Y,A10) 

I? (NU’IPRT.EQ. 2) WRITE (1,120) AKAME.HSAD.HEADl ,HEAD2 
FORMA T (" » ,A 1C , 3 X, A 1 0 , 83X, Ai 0 , 27 , A * 0) 

DM A X“ DAT A I 1 ) 

DMI 7= DATA 1 
DO 5C M-2 , MAX M 
DAT A. 1=0. AT. A (,’i ) 

IF ( DATAM.LT. DM AX) GO TO HO 
DM AX -DAT AH 

I ? ( DATA M * G T. DMIK ) GO TO 50 
DM I ,N=D A T A M 
CON . IM’J ! 3 
W ID?H=D 1.AX- DM IN 
DO 80 M=1,:i\XTJM 

K-1 . 0+ ( (DATA (M) -D w IN) /WIDTH) * TTY HUM 
tp f tf . _ 1 1 3Q 7 5 


I* (K. CQ. 1 ) 
IZ0D1=K-1 ' 

DO 


70 N=1 ,13 0^1 
C S ? 0 ?. H ( N) = ' r -" 

csrop .3 (K) ="+" 

IZ002-K+ 1 

DO 77 J=IZ 002 ,81 

C3T0RE (J) =" 

■J 5TRS 


I F ( !I r J 1PRT , EQ. 0) WRITE (1 , 9 r ) TI KES ( M) , DATA (M) , 
FORMAT (" 1?S1 A .4,1X, 1 PE1 1 . 3 , 7 X,.8C A 1 ) 

IF ( MUMPRT.FQ. 1) HP ITS (1 ,92) TI LIES (M) , DAT A (U) , 


DA? A 1 (M) 


(C STORE ( L) , L— 1 , 
(CSTORfi (L) , L=1 , 


I F ( NO U P ’R ? ♦ 3 Q • 2) WRITE (1 ,92) TIMES (M) , DATA (M) , (CSTORE(L) ,L=1 , 
DA r A 1 (M) , DA T A 2 ( M) 

FORMAT (" ", 1?E10. 4,1X, 1PE1 1 .3,1X,80M , 1X,1PBH. 3,1X, 1 PE1 1 
CONTINUE 
RETURN 
END 


M 

8 °) 

3 n ) 

• V 
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IV. DISCRETE DESCRIBING FUNCTION OF THE COMBINED 
WIRE-CABLE TORQUE AND FLEX-PIVOT NONLINEARITY 

In the last three chapters the continuous-data IPS system with the 
combined nonlinearity of the wire-cable torque and flex pivot is studied. 

The describing function analysis is applied to the continuous-data system. 

In the present chapter the digital IPS system with the combined 
nonlinearity is considered. The block diagram of the system is shown in 
Fig. *i-l. The combined nonlinearity is represented by the block N. The 
digital system is characterized by having the sample-and-hol d units which 
have the notation S/H. In order to apply the discrete describing function 
technique, a sample-and-hold .uni t is inserted in the path of the nonlinearity 
N. 

The objective of the investigation is to study whether self-sustained 
oscillations will occur in the digital IPS system. If self-sustained 
oscillations do occur, as they most likely will in a digital nonlinear 
system, what are the amplitudes and frequencies of these oscillations ? 

*1.1 Discrete Describing Function of the Combined Nonlinearity 

The discrete model of the combined nonlinearity of the wire-cable 
torque and the flex-pivot characteristics is shown in Fig. *1-2. The zero- 
order hold at the input of the combined nonlinearity is deleted, since 
there is already a zero-order hold at the output of N. 

The wire-cable torque T is modelled as shown in Fig. *1-2 and is 
functionally represented by 

T wc =H WT SGN( °i> +K V,T 0 i 
or 

T «< 0 i> = h wt + \n Q i ®i 2. 0 (,, - 2) 
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Figure 4-2. Discrete model of the combined nonlinearity 
of the wire-cable and flex-pivot torques. 


XX - - H wr + \n e i 


e. < o 

t — 


(4-3) 


where H lIT is in N-m, in N-m/rad. 0. is in rad, and T (0.) in N-m. 

WT wT i wc i 

It has been established that the Dahl solid rolling friction charac- 
teristics can be approximated by the nonlinear relation, 


where 


For 


dT Fp( 0 i) 

d©. 


= Y(T fp| - T fpq ) 


i = posi ti ve number 
Y = positive constant 

T FPI = T FP SGN i ) 

Tp P Q = saturation level of Tp p 

2, Eq. (4-4) is integrated to give 

© + c = "y" (t^~ t ) — - 

yw FP - FPO' 

i 2 Y(T fp + t fpo ) 


0 . > 0 
i — 


0 . < 0 


(4-4) 


(4-5) 


(4-6) 




where Cj and C 2 are constants of integration, and 


T + = T 
FP ' FP 

0. >_ 0 

(4-7) 

T“ = T ' 

'FP FP 

0. < ° 

(4-8) 


The constants of integration are determined at the initial point where 


0j . = initial value of 0. 
Tpp. = initial value of T pp 


Then, 


C 1 ~®ji " v(T + -~T J 
T ' 1 FP i FPO' 


C 2 "®i i ~ y (T+ . + T J 
11 YU FPi FPO' 


0 . > 0 


0j < 0 


(4-9) 
(4- 1 0) 


The describing function analysis depends on the assumption that the- 
input to the nonlinearity is a sine wave. Let 0.(t) be described by the 
cosinusoidal function, 



0. (t) = Acosoot 


(4-11) 

Then, 

• 

0. (t) = -Aws ino)t 

1 


(4-12) 

Thus, 

0 

11 

1 

> 

0. > 0 
1 — 

(4-13) 


0. . - A 

1 1 

0. < 0 
t — 

(4-14) 

The constant 

of integration in Eqs. (4-9) 

and (4-1 0 ) become 



^1 " A y(T + - T ) 
T v 1 FP i TP0 ; 


(4-15) 


r a x 1 


(4-16) 


A y(T + - T ) 
Y ' FPi TP0 ; 



Substitution of Eqs. (4-11) and (4-15) in Eq. (4-5) and simplifying, the 


solution of Tjjtp. is written as 



30 


T+ (t) = T 

*FP V ' 'FPO a 


R a 

R - _ Y + j(l - coscot) 


(1 - coscot) + 


1 


2 N ' ' R - l 

which is val id for 0. > 0 or (2k + 1 ) tt <_ cot <_ 2(k + 1 )tt, k 
where 


(4-17) 


0 > i i 2 , . . . i 


a = 2yAT 


FPO 


(4-18) 


R ‘ 


z 

a + 1 


FPi 

FPO 


(4-19) 


Similarly, for 0. £ 0, using Eqs. (4-11) and (4— 1 6) in Eq. (4-6), we have 


~ §0 - coscot) 


T - = t 

FP 'FPO a„ , 1 

-(1 - coscot) + 


(4-20) 


2 \ ■ *•«■*>*>«./ ■ R + | 
which is valid for 2kir £ cot £ (2k + 1 ) tt , k = 0, 1, 2, . . . 

The relations for T+p and Tpp in Eqs. (4-17) and (4-20) together with 
those of T^^.(0.) in Eqs. (4-2) and (4-3) are used for the derivation of the 
discrete describing function of the combined nonlinearity. As shown in Fig. 
4-2, the total torque disturbance due to the combined nonlinearity is given 
by 

T+ _ T+ + OJ+ 

n 'WC FP 


r- v + §•(!- coscot) 

= H WT + K wi Ac0SWt + 3T fpo a7 ( 1 

2-0 - coscot) + p — 


©j > 0 (4-21) 


T n T WC + 3T FP 

d ~ + 1 -7(1- coscot) 

” h wt + K w T Acosajt + 3T fpo ft; ; — 

7(1 - coscot) + 7 — - j - 

For the discrete model, we let 


0. £ 0 (4-22) 


0. (t) = Acos (cot + <j)) 


(4-23) 


where <p denotes the phase in radians. The z-transform of 0. (t) is 



(4-24) 


0.(z) = l Acos + <f>) z" 


k=0 


N = 


2,3 Or, in closed form, 

/ v Azf(z - cos27r/N)cos4) - s in27r/N*s info] 
! ' z - 2zcos2tt/N + 1 


(4-25) 


The z-transform of T p (t) is denoted by T p (z). Then, the discrete describing 


function of the combined nonlinearity is defined as 

T (z) 

N(z > - 67W 


(l|-26) 


The discrete describing function (DDF) for N = 2 is derived separately 
in the following section. 

4.2 The DDF of the Combined Nonlinearity For N = 2 

Let T p (kT) denote the value of T p (t) at t = kT. For N = 2, the signal 
T*(t) is a periodic function with a period of 2T. The z-transform of T*(t) 
is written as 

T n (z) = T„ ( 0) ( 1 + z -2 + z _2} + . . . ) + T n (T)(z _1 + z" 3 + . . . ) 


T p (0) tl + T n (T)z 


(4-27) 


z 2 - 1 

For the cosinusoidal input of Eq. (4-23), the corresponding expression 
for Tp p ( t) and T“ p ( t ) are 


T F P ft) - T 


~"T + y(l ~ COs(0Jt + (j))) 


FPO 1 


R - l 


+ |-(l - cos(o)t + 4>) ) 


0 . > 0 


t fpo 


R 


p - |{l - cos (tot •+ 4>) ) 


-(l - cos (cot + <}))) 


0 . < 0 
i — 


(4-28) 


(4-29) 


R + 1 ’ 2 

respectively. For t = kT, the last two equations become 
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T£p(k T ) = t fpo 


R - 1 


+ fO - cos(^p- + <(>)) 


+ &! -cos(^*i + <j>)) 


R - 1 


0 . > 0 

i — 


T FP< kT) = T FPO 


Tf- |0 ~ cos <ir + 

1 + §-(l - cos(^- + <j>)) 


R + 1 ’ 2 ^ N 

For N = 2, Eq. (4-25) is simplified to 


0 . < 0 

i — 


e.U) 


Substitution of Eqs. (if“32) and (4-27) into Eq. (4—26) , we have 

T (z) T n (0)z + T (T) 

N(z) ~ 07(17 " A(z - l)cos<{)' 


Also, for N = 2, z = -1; the last equation is simplified to 


N(z) 

For N = 2, 


T n (0) - V T) 

2Acos<{> 


T n (°) - T-(0) 


= T + (0) 


0 £ <J) < IT (0. £ 0) 
7 T £ <j) < 2 tt ( 0 . 5 0 ) 


T n (f) - T+(T) 


- ^ 


0 £ <p < n (0. £ 0) 


7 T <_ <f> < 2 TT (0. £ 0 ) 


Then, Eq. (4-34) becomes 

T"(0) - T+(T) 


N(z) = 


N(z) = 


2Aco5$ 

T+(0) - r(T) 


2Acostj) 


0 £ <p < IT 
7T < <j) <2lT 


For stability analysis, it is convenient to define 


(4-30) 

(4-31) 

(4-32) 

(4-33) 

(4-34) 


(4-35) 

(4-36) 
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F(z) = - 


N(z) 


Using Eqs. (4-2) , (4-3), (4-30) and (4-31), we have 


T n (0) = H WI + \JT AC0S * + 3T FPO “1 


R — _ • - j - + fo - cos$) 


+ ~( 1 - COS(j)) 


T n< 0) ■ " H WT - “w/' 05 * + 3T FP(T1 ' 2 


(1 - cos<j>) 


R , + 2"( ] “ cos ^) 


T n (T) = H WT + K WT Ac0s((f + 7r) + 3T FP0 — T 


■“ r + f (l - cos (tt + (j>) 


R • _ y + |(l - cos (tt + (j)) 


pT — jf + f(l + cos cj>) 

= h wt - *«*“♦ + 3T fpo 


R V - f ~ f(l + COS^) 

T " <T) = ‘"WT ' KwtAC ° S ‘* + 3TfP0 ' + f(, + cosTj" 

R + 1 2 


Thus , for 0 £ <j) < tt, 


T c (0) 


TJ(F) 


2K^.Acos4) 


2H wt - 3T fpo 


ITT- r - 2° ~ COSl)l) 

irrr + 1° ' cos ^ 


~-j- + “(I + cosij)) 


FTT + I (1 ' cos,|>) 


For ir j< cj) < 2tt, 

1^(0) - T”(T) = 21 ^^ 054 , + 2H wt + 3T fpo 


' 

R 


- COStj)) 

K 

- 

1 2 U 


1 

__ . a /1 

- cos<jO 

X 

- 

1 + 2 U 


TJTj* * jO + cos 40 

1 3 

— \ + ‘2^ 1 + cosc ^ 


(4-37) 


(4-38) 


(4-39) 


(4-40) 


(4-41) 


(4-42) 


(4-43) 


34 


4.3 Properties of F(z) - -1/N(z) for N = 2 as A 4 0 and A -► 00 

The properties of F(z) for N = 2 as A 4 0 and A 4 00 are now investigated. 
These properties will be useful in the determination of the critical regions 
of F(z) for stability studies. 


Thzotiw 4-1 . 


For N - 2, 

lim F(z) = 1 i m -1/N(z) = -1/Ky.j. for all <f> 

A-*» A- 400 


Proof: From Eq. (4-19), 


(4-44) 



R 


-1 + 

/a 2 + 1 


R 


1 

-1 - 

a +v^~ 

+ i 


R 


-1 + 

/^Tl 

- 

R 

+ 

1 

-1 + 

a +/a 2 

+ 1 


1 



a 


R 

“ 

1 

-1 - 

a +/a^~ 

+ 1 


1 



a 


R 

+ 

1 

- 1 + 

a +/a 2 

+ I 


,(4-45) 


(4-46) 

(4-47) 

(4-48) 


and a = 2yAT 


FP0‘ 


Let F + (z) = F(z) for 0 £ <p <n. Substituting Eq. (4-35) into Eq . (4-37) » 
we have 

F+ <*> - f-isr^r) 


Then, 


•2Acostj) 


lim F + ( 2 ) = lim p -fo r' - ~ T<TT) - 


(4-50) 


A 400 A 400 n 

Substituting Eqs. (4-39) and (4-40) into the last equation and using Eqs 


(4-45) through (4-48) , we get 
















Tneo/te/n 4-2. 


For N - 2, 

1 im F(z) = 1 i m 
A*0 A+0 


°/ 180 ° * iff; 


for all <j) 


Proof: For 0 < <}> < ir, 


1 im F + (z) « 
A-*0 


1 4Va*+ 1 


-1 + a +/a^+ 


1 + a +/a ,i + l 


+ 2K w _Acos<|> - 2H 


FPO 1 


TTT ~ 2 (1 ~ cos ^ _ TTT + 1° + COS(|)) 
TTT + 1° ' cosi>) ^TTT + f (l + cos 


, •. -2Acos4> 

1 im — rrj — 21 

A->0 WT 


= 0/180° 


for 0 <_ <f> < tt/2 


for 7r/2 < <J) <_ it 


Similarly, for tt <_ <p < 2n, 


1 im F“ (z) = 


+/a z + 


( 1 + COS(()) 
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1 im F“(z) 
A+0 ' • 


= 1 im 

"" A+0 


? t fpo 


-2Acos(|> 


■1 + a 


I 


+ “(i 


COS(J>) 


1 + a 2 


(1 + cos<j>) 


1 


[ 1 + a + 2 (1 


- cos<£) 


TT~a + 2 ^ + COS( ^ 


2K wt Acos(})+2H wt 


= 1 im 


-2Acos<t> 


A+0 2K wt Aco S {{) + 2H 


WT 


= 0/180° + tan" 1 j~f|y = 0/1 80° 3 tt/2 < <{> <2tt 


= 0/0 


7T <J) <3it/2 Q.E.D, 


4.4 The DDF of the Combined Nonlinearity For N > 3 


A general relation can be derived for the combined nonlinearity for 
N > 3. 

In general, the z-transform of the output of the combined nonlinearity 
may be wr i tten as 


CO N - 1 

T U) = I l T (kT)z 

M _a r a 1 * 


-k-mN 


N-k 


m=0 k=0 
N-l 
l T (kT)z 
k=0 n 

z N - 1 


The discrete describing function then becomes 


N-l 


T.(z) 


l T (kT)z' 


N-k 




k=0 


(z N - 1) l (Acos pp- + <p) 2 k 
k=0 N 


The last equation is simplified to 


N(z) = 


N-l 

I T (kT)z 
k=0 n 


N-k- 1 


A V( cos *)z* 

k=0 N 


(N > 3) 


(4-52) 


(4-53) 


(4-54) 


Or, alternately, 



(4-55) 


N-l 


\ • k— 0 

N(z) s n-7 


l T (kT) Z "- K (z 2 - 2zcos^ + I) 


][ (z - e^ 27rk/ ^ N )Az((z - cos-^-) cos<j) - sin^sincf)) 
k=0 


Or, 


N (z) = 


N-l 

I T (kT)z 
k=0 n 


N 

N-k-1 


N-2 


A(z - 1) l (z - e^ 2 ^ k ^)((z - cos^-)cos(j) - sin~sirx})) 


k=2 


For N - 3, z = ej 27r/3 ; 


(4-56) 


T (0)z 2 + T (T) z + T (2T) 

- n _ . t n - n _ 

' ' A(z - l) ((z + 0.5)cos<|) - 0 . 866s 1 n<f>J 


(4-57) 


Similar expressions can be obtained for N = 4, 5, • • . with z identified 
j2ir/N 


with e 

In general , 


T n (kT) - T-(kT) 


= Tj(kT) 


0 < ^ + * < it 

2TTk 

7T < ~ + $ < 27T 


(4-58) 


4.5 Asymptotic Properties of -1/N(z) for N > 3 as A Approaches Infinity 

The' fol lowing properties of -1/N(z) are found for N 3- The proofs 
of these properties can be obtained the same way as those for N = 2 by 
replacing 4> by <p + 27rk/N. 


(a) 

1 im 

(-1/N(z)) = - -L- 

for 

all 

<P 

(4-59) 


A 

\n 





(b) 

1 im 

|-1/N(z) | = J— 

for 

all 


(4-60) 


A-* 30 

Vr 




(c) 

l im 

(Arg(-1/N(z)) = 180° 

for 

all 

<P 

- (4-61) 


A-x» 





Asymptot 

:ic Properties of -1/N(z) 

for N"> 3 

as 

A Approaches 

Zero 


The following properties of -l/N(z) for N ^3 are obtained as A 



39 


approaches zero. 

(a) For N ^ 3> and for all <f>, 

lim | 1 /N (z) | = 0 (4-62) 

A+0 

(b) For N >_ 3, (N = odd integers) 

lim Arg(-1/N(z)) = - + - 2k - u + cp (4-63) 

A+0 Z 

for kTr/N < <p < (k + 1 ) tt/N , k = 0, 1,2, . . . , (2N - 1) . 

(c) For N :> 4, (N = even integers) 

lim Arg (-1/N(z)) = - - 2 - + -ylv - - tt + <J> (4-64) 

A+0 

for 2kfr/N £ cj) < 2(k + 1 ) tt/N . 

In view of the properties of ~1/N(z) listed above, the following theorems 
are generated. 

Tktotim 4-3. 

For even integral N >_ 3, the magnitude and phase of -1/N(z) repeat for 
every cj) = 2tt/N radians. 

Theosiejn 4-4. 

For odd N (N >_ 3) > the magnitude and phase of — 1 /N ( z) repeat for every 
= tt/N radians. 

4.7 Discrete Describing Function Plots of the Combined Wire-Cable and 

Flex-Pivot Nonlinearity - The Critical Regions 

The discrete describing function, N(z), for the combined wire-cable 
and flex-pivot nonlinearity is derived- in the preceding sections. The plots 
of F(z) = -1/N(z) together with the plot of G (z) , which is the linear 
transfer function that N(z) sees, in the frequency domain allow the study of 
the condition of self-sustained oscillations of the digital IPS system. 
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Computer programs for the evaluation of the -1/N(z) for N = 2 and 
N ^ 3 have been prepared. The listings of these programs are given in Tables 
4-1 and 4-2, respectively. 

For N = 2, the expression for F(z) = -1/N(z) is given in Eqs. (4~35) , 
(4-36) and (4-37)* Figure 4-3 shows the F(z) plot for N = 2 in the gain- 
phase coordinate's with 0 A < °° and all values of <p. 

The following set of parameters are used for the nonlinear elements: 



0.00225 

9.2444xl0 4 


V - 100 

h wt = 1 


In view of Theorem 4-1, Eq. (4-44), the most important parameter among those 
listed above is , since when A approaches infinity the magnitude of F(z) 
approaches l/K^. However, as shown in Fig. 4-3, the F(z) plot stays on the 
-180° and -360° axes for N = 2. 

The discrete describing function N(z) for N ^ 3 is given by Eq. (4-54) 
or Eq. (4-56). Figure 4-4 shows the gain-phase plot of F(z) = -1/N(z) for 
N = 3- The curves for several values of tj) between 0° and 60° are plotted to 
illustrate the effect of varying the phase of the input signal to the non- 
linearity. It should be noted that for N = 3, Theorem 4-3 states that the 
values of F(z) repeat every 60 degrees starting from <j> = 0°. As the magnitude 
of the input signal, A, approaches infinity, the magnitude of F(z) becomes 
1/K^.j. which is -40 db in this case, since is 100. On the other hand, 
as A approaches zero, the magnitude of F(z) becomes zero or - 00 db, and the 
phase of F(z) is bounded by -300° and -24 0° for all values of ( J> . When the 
value of K^-j. is varied, the curves of F(z) will shift up or down according to 
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Table 4-1. Computer program, for the computation of the discrete describing 
function of the combine flex-pivot and wire-cable nonlinearity 
of the IPS. N = 2. 

TT CALCULATION FOR ~1/N(Z> FOR COMBINED NONLINEARITY y N-2 

COMPLEX GN1 vGN2 

REAL*8 P < 20 ) y RAD y TO y GAMMA t AST ART y A y AA y R y CSS y CSP r CSC 
~ REftL*8 IT A 2 y Rl y R2 y R3 y R4 r TOP y TOM y TTH y TIP y T1 y 12 

PI -3 ♦ 14159D0 
RAD- I 80 ♦ DO/PI 

“ — — CSP-T7TR5 : 


T0~0 ♦ 00225D0 
KWT-IOG « DO 

TIWT-1 .DO 

GAMMA-- 9 ♦ 2444D4 
ASTART-1 .OD-TO 


RKPBODDQi; 


NP-U 

ND-15 

WRITE < 5 y 100 ) 


■ write 'rsrrorr 

DO 1 J— 1 y ND 
DO 1 I-lyNP 

css«i*Do-csp ‘ 

CSC-1 ♦ DOtCSP 

A-ASTART*DFLOAT < I > * < 10 ♦ DO** < J~ 1 > ) 


Aft- 2TIJTO3 A M M ft TO 

R- ( -1 ♦ D0/AA > +DSQRT ( < A A* A Ail ♦ DO > / < A A* A A ) > 
AA2-AA/2 ♦ DO 


RT-I . DO/ ( K'-l ♦ DO ) 
R2-R*R1 

R3-1 ♦ DO/ ( Ril ♦ DO ) 


R4-R*R3 

TOP- < R2+ AA2*CSS ) / < R 1 iAA2*C8S ) 
TOM- ( R4“AA2*CSS > / ( R3iAA2*CS3 ) 


f I M- ( R4~AA2*CSC ) / ( R3iAA2*U8in ! 

TTP- < R2i AA2*CSC > / < R 1 i AA2*C8C ) 

Tl-3 . DO* TO* ( TOP “TIM ) +2 ♦ DO* ( HWTiKWT*AA*CSP ) 


T2-T;XiOTrO* ( TOM -TTP 5 +2 ♦ DO* H~iWTiKWT*AA*CSP > 
GN1— 2 ♦ D0*AA*CSP/T1 
GN2--2 ♦ D 0 * A A * C S P / T 2 


Gll-REALTGNTT - 
G12-AIMAG < GN1 > 

e w\ <1 .... i”. r" a i / v i /*> \ 
i A-:. 1 — TK C- H L \ b ^ i 


G 2 2 - ft' I nft 'G ( G N 2 ) 
GMAG1-CABS < GN1 > 
GM AG 2-CABS ( GN2 > 


■BDB1-20 , *ALGG10 < GMAG1 ) 
GDB2-20 . *AL0G10 < GMAG2 ) 

G P H 1 - R A D * A T A N 2 < 012 y G1 1 > 


n ! ! n . .. P'l *■ T*V“i;"A-»V" _ A-AT."S / I'SW fVWrT 

IF ( GPH2 . GE .0. > GPH2-GPH2-360 . 
IF (GPH1 *GE«Ot > GPHl-GPHl“~3a0 ♦ 


'WRITE ( 3 r .1.02 7ft y GDB2'y GMAG2 y GPl-12 
CONTINUE 

FORMAT < * DISCRETE DESCRIBING FUNCTION FOR IPS ') 
FORMAT w' r SX y ' A ' y 1 1 a y 7 GDB2 ' y 10 a y ' GMAG ' y 9X y r PHASE ' / ) 
FORMAT ( 1P4E14 ♦ 0 > 

STOP 


1 

100 

tot 

102 



Table 4-2. Computer program for the computation of the discrete describing 
• function of the combined flex-pivot and wire-cable nonlinearity 
of the IPS. N GE. 3- 


DISCRETE DESCRIBING FUNCTION FDR COMBINED NONLINEARITY? H. 
REAL PHI *P< 15> -PI -RAD - GAMMA - ASTART -A ?AA »R *HWT >KWT >T0 
REAL PP $ RP - AAR - R1 - RS » R3 - R4 > A fl£ - TC j PH I K ? PH I D - p I K -TUN - TUP 
COMPLEX GV jTTSUM 

COMPLEX TSUMN - ZSUMN * ZSUM - T2 > TSUM ? 2 » THETA - GNN - GN 
PI =3. 141 59 


RAD -ISO ..•-'PI 
HUT=1 . 0 
KMT = 1 0 0 . 

GAMMR=9 .'2444E4 
AST ART = 1 .E-8 
ND=9 
HP-1 

TO=0 . DOES 5 
RR = 0 . 

N I =50 

AN=FL0AT CNI ) 

N1=NI-1 
N2=NI -£‘ 

P<1)=1 .0 
NP HI =36 00 
PP=2 .♦PI ■••'FLOAT ONI ) 

T.HETfi=CMPLX C RR j PP ) 

Z=CEXPC THETA ) 

ZSUMN=CMPLXC1 . - 0 . 0> 

TS’JMN=CMPLX i 0 . 0 - 0 . 0 > 

WRITE <5 - 1 00) 

MR I TEC 5 * 1 02 ;> AN 
WRITERS j 1 1 0) GAMMA 
MR I TEC 5 fill > HWT 
MR I TEC 5 , 1 1 £) KMT 
MR I TEC 5 j 1 1 3) TD 
MR I TEC 5 * 1 01 ) 

DO o I = U ? 1 44 j 36 
PHI-C2 . 0*P I ♦FLOAT C I > >X C FLOAT CNPHI > ) 
PH I D=PH I ♦RP.D 
MR I TEC 5 j 103) PH ID 
DO 1 J=1 ?ND 
DO 9 L= 1 -NP 
2SUMN=CMPLXC 1 .0*0.0) 

T S ! JMN=CMPLX C 0 . 0 j 0 . 0 ) 

a=astart ♦float c l > ♦ C 1 0 .♦♦C-J-l )) 

HA - 2 . OtGAMMA^A^TO 
A AS = A A- . 0 



Table k-2. (Continued) 


R = C-l . 0.--RR > +SQRT < < RRaRR+ 1 . 0 > C ARARR > > 

Rl = l .0- R-l . CO 
R2-RARI 

R3=l .0/CR+1 . 0 > 

R4=R*R3 
DO 2 K= 0 s N 1 

PIKX2 . OaP I AFLOAT OO VCFLOAT CNI > > 

PHIK=PHI+PIK 

I F C PH I K . GT . C 2 . UAP I > :> PH I K=PH I K-2 . 0*P I 
TWN=-HWT+KWTaR aCOS C PH I K > 

TUF-HUT +KUT a RaCOS C PH I K > 

IFCPHIK.LT .PI >GQ TO 6 

TC=3 . OATD a C R2+flR2-Rfl2ACDS < PH I K > ) •••' C R 1 + RR2-RR2ACDS C PH I K > > +TWP 
GO TO 7 

TC=3 . OaTOa C R4-RR2+RR2aC 0S C PH I K > >^<R3+fiR2-fifl2ACDS<PHIK> > +TUN 
T.SUf'1=TC aZa a H 1 -K :> 

T S U M N = T S U M N + T S U M 
CONTINUE 

IF CHI .LE .3 >60 TO 10 

DO 3 M=2 ? N2 

ZSUM=Z-ZaaM 

ZSUMN=ZSUMNaZSUM 

CONTINUE 

CONTINUE 

TZ= C Z-COS C PP > > aCOS C PHI > -3 1 N C PP > AS I N C PH I > 

6N=TSUMrV C RaTZa < Z- 1 . > aZSUMN > 

GNN=-1 . .•■"GM 

6V'=GNN 

H 1 —REAL C GV ’> 

G2=RIMAGCGV> 

G M A G — C R B S C G V ? 
bBB=£ U . ARLOb 1 U C bMRG ) 

GPHAsE— RAD aRTRNE CG2 >G1 ? 

IFC bPHASE . bE . U « .•'GF'HASE = GFHAS'E~360 . 

UR 1 1 E C 5 ? 1 04 R j bPHRSE j bDB * bMRG 
FORMAT C1PE24.15) 

CONTINUE 

CONTINUE 

CONTINUE 

FORMAT C 7X BE SCR I B I NG FUNCT I ON GF COMB I NEB NGNL I NEAR I TV •" > 
FORMAT C 9 ?X j R 1 OX j •" PHASE j 1 OX » BE ' j 1 OX j MRGN I TUBE ' > 

FORMAT C 5 X ? N = ,F4.i ;• 

FORMAT C ,5X > ' GAMMA = •" s 1FE12 ,4> 

F □ RM H T C 5X j H U T= ,- ‘ >F 3 . 2 > 

FORMAT C SX j •" KUT= •" j FS ,£} 

FORMAT C OX 9 •" TD= y » 1 PE 1 £ . £ > 

FORMAT C j 5X , PH I - •" ?F5.1 > 

FORMAT C 1 P5E 1 4 . 5 > 

STOP 

END 



±n 
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1/K^.j. in db; when the values of TppQ> Y> and H^. are varied, the shape of 
the curved portions of the plots in Fig. 4-4 will be changed. However, in 
general, the impact of the variation of K^. will be the greatest. 

The F(z) plots for N = 4, 5,6,8, 10, 20, and 50 are shown in Figs. 4-5 
through 4-11, respectively. For N = 4, the F(z) plot extends from “315° to 
-225° as (p varies. For N = 5, the span of the plot is -288° to - 252 °. 

For stability analysis, it is sufficient to consider only the bounds of 
the F(z) plot for a fixed N. Self-sustained oscillations characterized by 
the period T c = NT, where T is the sampling period, may occur if G e ^(z) 
intersects with any part of the F(z) plot. The region bounded by all .the 
F(z) curves for a given N is defined as the critical region . The critical 
regions for the combined nonlinearity for N > 2 are the regions that are 
bounded by the F(z) curves for (p = 0° and tj) = 2tt/N for N = even and c|) = tt/N 
for N = odd. 

As N approaches infinity, the discrete describing function N(z) approaches 
the describing function N of the continuous-data nonlinearity, as shown in 
Fig. 4-1 i. It is observed that as N increases the width of the critical 
region becomes narrower. As N approaches infinity, the critical region of 
F(z) = -1/N(z) approaches the -1/N plot shown in Fig. 2-1. 
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.Figure ^-5. Discrete Describing function plots of combined flex-pivot 
and wire-cable nonlinearity for the IPS. N = k. 
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Figure *4-6. Discrete describing function plots of combined flex-pivot 
and wire-cable nonlinearity for the IPS. N = 5. 
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Figure k~7 . Discrete describing function plots of combined flex-pivot 
and wire-cable nonlinearity for the IPS. N = 6. 
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Figure 4 - 11 , Discrete describing function plots of combined flex-pivot 
and wire-cable nonlinearity for the IPS. N = 50 . 
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